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Oh' Discrete lattice simulations of a one-dimensional rf> 4 theory coupled to an external 

heat bath are being carried out. Great care is taken to remove the effects of lattice 
discreteness and finite size and to establish the correct correspondence between 
, simulations and the desired, finite-temperature continuum limit. 

1 Introduction 

m 

To be able to study numerically certain properties of cosmological first order 



phase transitions such as the nucleation rate of bubbles, we investigate first the 
effects that finiteness and discreteness of a lattice have on results derived from 
, simulations. For this basic study we limit ourselves to one spatial dimension 

Q\ ' and to a real scalar field 4>(x,t) subjected to a potential Vb(0) and an envi- 

ronmental temperature T. The dynamics obey a Langevin equation as given 
in Borrill and Gleisertl. This 2D-paper indicates that for a given tree-level po- 
■ , tential, results obtained numerically on the lattice can't straightforwardly be 

identified with analytical results. They employ a renormalization procedure to 
r"| , get rid of the lattice spacing dependence and to identify the correct continuum 

limit of the simulations. Similar problems in ID are treated in what follows. 



2 Method 



For classical field theories, the one-loop corrected effective potential is given 
by a momentum integral □, and evaluated to 

v 1L &) = v (<t>) + \£ gin [i + = km + - y / TO) ■ (i) 

The discretization 5x of the lattice and its finite size L introduce short and 
long momentum cutoffs fc m i„ = 2n/L and A = n/Sx. Therefore the simulation 
only sees V!l = Vo((f>) + (T/2) f k . . .dk. If one neglects the effect of fc m i„ 
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which is possible for sufficiently large L [J integration and expansion in powers 
of Vq /A 2 (possible for sufficiently large A) yields 

fe-W + I^-^f + AO^). (2, 

As is to be expected for a one dimensional system, the limit A — > oo exists and 
is well-behaved; there is no need for renormalization due to divergences. How- 
ever, the effective one-loop potential is lattice-spacing dependent through the 
explicit appearance of A, and so are the corresponding numerical simulations, 
as evidenced in the tree-level potential cases of Fig. 1 (left graphs). 

Similar to the renormalization procedure for 2D systems given by Borrill 
and Gleiser El we remove this dependence on Sx by adding counterterms to 
the tree- level potential V). In contrast to higher-dimensional systems, these 
counterterms are finite, namely Vct(</>) = (T '/ '47r)(V r " '((/>) / 'A) . Hence the lattice 
simulation works with the corrected potential 

TV"(rti\ 

VuttW = W) + -^^fo (3) 

and simulates the continuum limit to one loop, Vix,(<^) = Vq + (T/2) . . .dk = 
Vil(^) (where V^'att * s employed in the integrand), just as it should be. 



3 Application 

Since the numerical extraction of bubble nucleation rates is a contrived pro- 
cess the ideas of Section 2 are tested initially with the symmetric double well 
potential V (4>) = (A/4) (0 2 -<£§) . We compare simulations using Vb alone 
with those employing Vtatt(0) = Vq(4>) + 3XTSx<j) 2 /A-k 2 (eq. ||). One set of 
runs investigates the mean field value (4>) of the metastable equilibrium before 
the first kink-antikink pair occurs (4>(t) = {1/L) J (f>(x,t)dx). Another set of 
runs measures the kink-antikink pair density n p (proportional to the number 
of zeros of the low-pass filtered field). Fig. I shows the comparison for different 
lattice spacings Sx. Apart from a discrepancy for very coarse grids (Sx w I) 
the average field value is clearly lattice-spacing independent in the right panel 
(VLatt), m contrast to the use of Vq. The effect is even more striking in the case 
of rip where the addition of the finite counterterm removes any Sx dependence. 

^Lattice simulations only know one size parameter, the number of degrees of freedom N = 
L/Sx. With a given N it is always possible to choose L big enough for the effects of 5x to 
dominate over those of L. 
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In summary it was demonstrated that even in a field theory without di- 
vergences finite counterterms play a role. Their inclusion gets rid of the de- 
pendence of simulations on size and lattice spacing. This can be observed in 
the averaged field value and in the density of kink-antikinks. Further stud- 
ies of this renormalization procedure a identify the correct continuum limit of 
simulations, thus matching theory and numerical results. 
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